NASA Technical Memorandum 87292 


Influence of Fatigue Crack 
Wake Length and State of 
Stress on Crack Closure 

M 8 6-226 8-5 

Ujicla s 
059 1 1 


f2JASA-T£l-8729 2) IHJiSEHCE OF iAXIGGE CEACK 
8 MS LEUGSH A HO SIA1E Cl SIC L£S CE CBACK 
CLOSURE (21A3A) 21 o TO AO 2/11 AO 1 CSCL 11? 

G3/26 


^•ifack Telesman and Douglas M. Fisher 
Lewis Research Center 
Cleveland, Ohio 



Prepared for the 

International Symposium on Fatigue Crack Closure 
sponsored by the American Society for Testing and Materials 
Charleston, South Carolina, May 1-2, 1986 


NASA 


INFLUENCE OF FATIGUE CRACK WAKE LENGTH AND STATE OF STRESS ON CRACK CLOSURE 


Jack Telesman and Douglas M. Fisher 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

ABSTRACT 

The location of crack closure with respect to crack wake and specimen 
thickness under different loading conditions was determined. The rate of 
increase of K^ in the crack wake was found to be significantly higher for 
plasticity induced closure in comparison to roughness induced closure. 
Roughness induced closure was uniform throughout the thickness of the specimen 
while plasticity induced closure levels were 50 percent higher in the near 
surface region than in the mid-thickness. The influence of state of stress on 
low-high load interaction effects was also examined. Load interaction effects 
differed depending upon the state of stress and were explained in terms of 
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INTRODUCTION 


In recent years much attention has been focused on fatigue crack closure 

as introduced by Elber [1]. He suggested that crack surfaces can close before 

the minimum load Is reached due to the residual plasticity induced 

displacements In the crack wake. He further suggested that the fatigue crack 

propagation (FCP) rate is a function of the stress Intensity range when the 

crack tip is open. Thus, FCP is dependent on AK .. and not on AK . More 

eff app 

recently, a number of researchers [2-4] have shown that in addition to 
plasticity, crack closure can also be caused by the mismatch of fracture 
surfaces (roughness) or by a build up of the oxide layers on the fracture 
surface. 

The importance of crack closure in present day research on FCP can be 
seen by the application of this concept to the research of such diverse topics 
as variable amplitude loading [5], near-threshold behavior [6,7] and the 
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accelerated growth of short cracks [8,9]. However, many questions regarding 
the actual mechanisms of crack closure still persist. For Instance, the 
location and distribution of closure stresses In the crack wake Is not well 
understood. It has been suggested that crack closure Is a function of crack 
wake length [10]. Zalken and Ritchie [11] as well as others [12,13] have shown 
that In the near-threshold regime where roughness Induced closure Is Important, 
removal of crack wake reduces closure. However, the Influence of crack wake 
Is unknown at higher stress Intensities where plasticity Induced closure Is 
dominant. To further complicate the picture, surface effects due to localized 
plane stress conditions can also significantly contribute to crack closure. 
Thus, in order to better understand the different closure mechanisms, there is 
a need to determine how and where closure occurs. It is also Important to 
understand the mechanism under which predominantly roughness or oxide Induced 
closure changes to plasticity Induced closure. 

Crack closure may also be Important In understanding load sequence 
Interaction effects which are not well understood. Review of the literature 
reveals some controversy regarding the FCP behavior Immediately following a 
low-high loading sequence. VonEuw et al. [14] observed a transient region 
after the transition to a higher AK where initial FCP rates were higher than 
expected. However, others [15,16] concluded that no load interactions take 
place when this type of loading Is applied. A very Interesting finding 
relative to a low-high loading sequence behavior was documented recently by 
Schulte et al. [17]. They also noted a region of higher than expected FCP 
rates after the transition to a higher AK when FCP rates were measured at 
the specimen surface. However, when the FCP rates of the same specimen were 
evaluated by using a striatlon counting technique in a region away from the 
surface, no initial acceleration of FCP rates was detected. The above 
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observations suggest that plane stress and plane strain variations should be 
considered in the effort to explain loading sequence effects. 

A test program was developed to consider the effects of crack wake and 
specimen thickness on crack closure, as well as the influence of the stress 
state on low-high load interactions. 

EXPERIMENTAL PROCEDURE 

For the FCP studies, compact tension ( CT) , specimens oriented in a long 
transverse orientation were machined from a 3.2 mm thick sheet of 7075-T6 
aluminum with tensile properties as detailed in Table 1. The width of the CT 
specimens was 38.1 mm, the initial a/w ratio was 0.2, and the thickness was 
that of the as received sheet unless otherwise indicated. All tests were 
performed using a closed loop servohydraul 1 c fatigue machine. All testing was 
done at an R ratio of 0.1 and frequency of 20 Hz in an ambient air 
environment. Testing was computer controlled using a compliance technique to 
maintain either a constant stress intensity range or load range and to 
determine crack length. Occasional optical readings were made to confirm the 
compliance readings of the crack length. 

Constant load amplitude tests were performed on the 3.2 mm thick specimens 
to obtain baseline data. Two types of tests were conducted, wake removal and 
surface removal, to determine the effect of crack wake length and thickness on 
crack closure. Initial preparation of specimens for both test types was 
precracking under constant stress intensity ranges of 2.5, 8.8, and 17.6 MPa/m 
until an a/w of 0.4 was reached. For the crack wake removal experiments, 
increments of the fatigued crack wake were removed using a 0.5 mm thick 
jewelers saw, and closure loads determined. The closure loads were taken at 
the first deviation from linearity in the unloading portion of the load 
displacement curves. In the surface removal tests, specimen thickness was 
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reduced Incrementally by machining equal amounts from each surface. At each 
thickness closure loads were determined by the above mentioned procedure. 

The low-high load sequence tests to evaluate stress state influence were 
performed on as-received and 1.1 mm thicknesses. The low-high load sequence 
tests performed on the 3.2 mm thick specimens consisted of precracking 
duplicate specimens at a AK of 4.4 MPa/m from the initial a/w of 0.2 to 
0.35. The stress intensity range was then increased to a AK of 13.4 MPa/m. 
In addition, separate CT specimens were tested at a constant AK of 4.4 MPa/m 
and 13.4 Mpa/in. For each of these constant AK specimens load-displacement 
curves were obtained, from which closure loads and the corresponding K CL 
values were determined. For the closure reading tests the frequency was 
decreased to 0.05 Hz. A similar procedure was used for the 1.1 mm thick 
specimens, however, the low-high loading sequence consisted of a transition 
from a AK of 6.6 MPa/m to a AK of 16.0 MPa/m. Again closure loads for the 
thinner specimens were obtained by testing separate specimens individually 
fatigued at the above two stress intensity ranges. The entire program test 
matrix is shown in Table 2. 

To determine the FCP rates as a function of crack length from the 
transition to a higher AK, f ractographic analysis of the striation spacings 

j 

was performed using a scanning electron microscope. All fractographs were 
taken using a 0° tilt to limit possible distortions of striation spacings used 
to calculate FCP rates. 

RESULTS AND DISCUSSION 
Baseline FCP Data 

The baseline fatigue crack propagation curve obtained for the 7075-T6 
alloy is shown in Fig. 1. An equation fitting the straight line portion of 
the curve (i.e., crack growth data for AK > 8MPa./m) is: 
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( 1 ) 


da , 
dN = 


,,, , n 10 Al ,3.07 

.236 x 10 A K 


where da/dN is the crack growth expressed in m/cycle and a K is in MPa/m. 

Removal of Crack Wake 

The data showing crack closure stress intensity^ as a function of crack 
wake length at aK of 2.5, 8.8, and 17.6 MPaTm is shown in Fig. 2. These 
results show a significant difference in behavior. The test conducted at the 
near-threshold stress intensity of 2.5 MPa/m, showed K^ decreasing when 
the remaining crack wake reached approximately 3 mm. The decreased 

gradually with further removal of crack wake. At the shortest crack wake for 
which valid measurements were obtained of 1.4 mm, K was measured to be 
0.75 MPa/m which was approximately 60 percent of the initial full crack wake 
value. These results agree with trends observed by others in both aluminum 
[11] and steel [13] where it was observed that in the near-threshold region 
K started decreasing when the length of the crack wake is reduced to 2 to 

l- 

4 mm. 

In contrast to the tests performed at AK of 2.5MPa/m, no change in 
K CL * he remova l of f he crack wake was detected for tests performed at 

AK of 8.8 and 17.6 MPa/m even though the shortest remaining crack wakes for 
which data was obtained were respectively 0.33 and 0.68 mm. Since there can 
be no closure without crack wake, the results indicate that closure load 
Increases with respect to crack wake much more rapidly at aK of 8.8 and 
17.6 MPaYm in comparison to a AK of 2.5 MPa/m. The divergence in behavior 
suggests that different closure mechanisms control closure at near-threshold 
stress intensities than those controlling closure at higher stress intensities. 


^For all conditions K 0 was found to be approximately equal to Kql 
thus AK e ff = K max - Kcl - 
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In the wear experiments, which were conducted to examine the effects of 

magnetic-head wear on magnetic signals, a newly recorded tape (11 m long) was 

first run and the standard signal was read back. An unrecorded tape 

(230 m long) was repeatedly run over a desired head sliding distance. During 

the sliding process the 11-m-long recorded tape, which was used to read back 

signals G at various stages of the experiment, was replaced several times, 
n 

Electron Diffraction and Depth Profiling 
The sliding surfaces of the magnetic heads were depth profiled to examine 
the crystalline states as a function of depth. The etching was done with 
hydrochloric acid at 50±1 °C for a predetermined etching time. The surfaces 
were examined by reflection electron diffraction in a transmission electron 
microscope. The acceleration voltage was 75 or 100 kV. The depth In terms of 
the etching time was obtained by using an optical interference microscope. 

RESULTS AND DISCUSSION 
Crystalline State of Lapped Surface 
Sliding a lapping tape on a Mn-Zn ferrite magnetic head generated 
abrasion and developed a deformed layer on the ferrite surface. Figure 3 
presents reflection diffraction patterns obtained from the abraded surface and 
the etched surfaces. The arcs In the electron diffraction pattern of the 
abraded surface are partial Debye-Scherrer rings that contain enlarged, 
streaked spots [Fig. 3(a)]. They Indicate that a highly textured, surflclal 
layer formed on the single-crystal magnetic head during the sliding of the 
lapping tape. 

The surface etched to a depth of 0.3 ym from the abraded surface shows an 
enlarged streak-spot pattern [Fig. 3(b)]. The streaking Indicates a great 
amount of plastic deformation: that Is, a highly strained single-crystal 

structure. A large number of line defects can cause streaking in a 
single-crystal diffraction pattern. 
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The surface etched to a depth of 0.6 pm shows Klkuchi lines (pairs of 
black and white parallel lines), which are an Indication of the bulk 
single-crystal structure of a ferrite head containing no extrinsic mechanical 
stress [Fig. 3(c)]. Thus the total thickness of the deformed layers on the 
Mn-Zn ferrite head presented In Fig. 3 was less than 0.6 pm. 

When lapping tape was again slid against the etched surface typified by 
Fig. 3(c), the polishing process reintroduced local surface stresses and 
redeveloped a deformed layer on the ferrite similar to that on the surface 
presented in Fig. 3(a). 

Effect of Crystal Structure on Magnetic Signals 
To determine how wear of the lapped magnetic head affects its crystalline 
state and magnetic signals, we ran the head against magnetic tapes (recorded 
and unrecorded) at an initial tape tension of 0.2 N and three head 
displacements for a sliding distance of 40 km. The results are presented In 
Fig. 4(a). The readback signal amplitude (sensitivity) decreased rapidly 
during the Initial 20 km of tape sliding distance. It remained constant or 
decreased slightly over the next 20 km. 

Typical electron diffraction patterns taken from the lapped surfaces of 
the magnetic head before the sliding experiment and from the wear surfaces 
after 40 km of sliding are presented In Fig. 4(b) and (c). These diffraction 
patterns show that sliding action changed the highly textured, polycrystall Ine 
structure of the lapped head to a nearly amorphous structure. When a lapping 
tape was then slid against the wear surface of the lapped magnetic head, the 
nearly amorphous structure was removed and a highly textured, polycrystalline 
structure was reformed In the magnetic-head surficlal layer. The lapping also 
restored the readback signal amplitude to the initial standard level G^ . 

In other words, lapping the wear surface of the magnetic head restored 
extrinsic magnetic properties. 
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with the trends observed In the crack wake removal tests which also indicated 
the rapid growth of closure for a plasticity dominant closure mechanism. 

Transition From Roughness to Plasticity Induced Closure 
The experimental results obtained in this study shed some light on the 
conditions under which the transition from roughness to plasticity dominated 
closure mechanism occurs. It was shown that even under a predominant plane 
strain condition encountered in the specimens tested at AK of 8.8 MPa/m, 
plasticity induced closure plays an Important role. Thus, the calculated 
overall state of stress of the specimen is not adequate for the determination 
of the type of predominant closure 
It can be postulated that 


mechanism. 


= K 


CLR 


+ K 


CLP 


( 2 ) 


where K is the closure stress intensity due to roughness Induced closure 
CLR 

and K CLp is the plasticity Induced closure. is dominated by the area 

in the crack wake where the largest mismatch of the crack faces occur. It was 

shown that for the specimen tested at A K of 2.5 MPaTm, the closure was 

uniform throughout the specimen thickness Indicating that K CLp was 

considerably smaller than K resulting in roughness dominated closure 

L LK 

mechanism. Due to the considerably larger plane stress plastic zone sizes of 
the 8.8 and 17.6 MPa/m specimens (Table 3), the near surface displacements at 
these stress intensities were also considerably larger in comparison to the 
test performed at the near-threshold stress intensity. This resulted in 
confining most of the closure near the surface and promoted plasticity induced 
closure. It is clear that the transition from predominantly roughness Induced 
to plasticity Induced closure occurs when K becomes a significant 

LLi 

component of the total measured K^. More research is required to obtain a 
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more quantitative description of the transition by measuring both the 
thru-thickness as well as near surface displacements at various applied stress 
Intensities . 

Low-High Loading 

While some researchers [15,16] have not found any load Interaction effects 

when a low-high loading sequence is applied, others [14] found Initially higher 

than expected FCP rates after transition to a higher AK. Hertzberg [20] 

proposed that these initially higher than expected FCP rates can be explained 

by the influence of the transitional behavior of crack closure on AK ... He 

ef f 

argued that initially after the transition to a higher AK, the material still 

experiences low closure loads and hence a greater AK^. Only when the new 

larger deformed zone created by the higher loads begins to interfere in the 

wake of the crack front, does the closure level begin to rise, decreasing the 

AK rr until a new equilibrium value of K_, and AK .... Is reached. This 
eff CL eff 

process is shown schematically in Fig. 5. 

In order to evaluate Hertzberg's hypothesis, baseline FCP rates as a 
function of AK eff were determined and are shown in Fig. 6. 

A Paris type equation fitting this data is: 

jjj = 1 . 445x1 0 _1 ° AKj f ^ (3) 

i 

where da/dN is in m/cycle and AK .... is in MPa/m. The least square 

err 

correlation coefficient for this equation Is 0.997. 

The results of the low-high loading sequence tests are shown in Figs. 7 
and 8. As shown in these figures, there exists a substantial difference in 
the behavior between the 3.2 mm and the 1.1 mm thick specimens. The thick 
specimens exhibited no noticeable change In FCP rates as a function of fatigue 
crack distance after the transition to a higher AK as shown in Fig. 7. 

However, for the thin specimens (Fig. 8) immediately following the transition, 
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the FCP rates were initially high and progressively decreased till a plateau 
was reached. For these specimens, the FCP rates Immediately after the 
transition were approximately twice as high as the rates In the plateau region. 
The above results were obtained by measuring striation spacings. Typical 
fractographs used for this purpose are shown In Figs. 9 and 10. 

The validity of using the striation counting technique to determine FCP 
rates was evaluated by comparing the growth rates obtained by this method to 
rates measured by the compliance technique. Striation measurements 
consistently showed slightly lower growth rates (up to 25 percent) possibly 
due to a small tilt of the fracture features with the respect to the screen 
projection. Some scatter within the striation measurements was also present. 
The scatter band was approximately ±10 percent of the mean value. This degree 
of scatter is reasonable when it is considered that a striation counting 
technique measures localized growth rates. These rates can be influenced by 
such factors as grain orientation and the presence of inhomogeneities such as 
inclusions . 

The question arises why there is a difference in behavior between the 
thick and thin specimens? The difference in observed behavior may be 
explained by evaluating the state of stress of each specimen and its effect 
on K . The applied stress intensities and specimen thicknesses were chosen 
to evaluate both plane strain and plane stress effects. The 3.2 mm thick 
specimen still met the ASTM E399 thickness criteria for a predominant plane 
strain condition even after the transition to a higher AK. However the 1.1 mm 
thick specimen was in a substantially plane stress state after transition. 

These differences in the state of stress had a considerable effect on K 

w L 

as is described below. 

Equilibrium K values for all applied stress intensities used in this 

U L 

part of the program are shown in Table 4. For the 1.1 mm thick specimens, the 
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K cl was measured to be 2.3 MPaVm at the lower stress Intensity and 5.3 MPa/m 
at the higher stress intensity. This substantial difference in K CL values 
is not surprising considering that the lower stress intensity resulted in a 
predominant plane strain state while the higher stress intensity resulted in 
existence of a substantial amount of plane stress. If it is assumed that 
Hertzberg's transition hypothesis is correct (Fig. 5), then the accelerated 
FCP rates immediately after transition to a higher AK can be predicted by 
the use of the da/dn-AK^ relationship (Eq. (3)). By substituting the 
appropriate AK e ff values into this equation, a predicted FCP rate of 
9.5xl0~ 7 m/cycle is obtained. The FCP rate obtained from striation spacing 
measurements, shown in Figs. 7 and 8, is lxlO -& m/cycle. The crack closure 
hypothesis (Fig. 5) also predicts that after a certain amount of crack growth, 
an equilibrium FCP rate should be reached corresponding to the baseline data. 
The predicted FCP rate at equilibrium using Eq. (3) is 4.8xl0~ 7 m/cycle. 

The actual FCP rate as measured from striation spacings is 5xlCf 7 m/cycle. 

Thus the crack closure hypothesis predicts the observed behavior very well in 
the case of the thin specimens. 

The acceleration of the FCP rates for the thin specimen, after the 
transition to a higher AK lasted for 0.2 mm of crack growth. It is 
interesting to note that the calculated plane stress plastic zone radius at 
the higher AK is 0.18 mm and is thus very close in size to the fatigue 
crack wake length over which accelerated FCP behavior was observed. 

Assuming that the initial acceleration of FCP rates after the transition 
to a higher AK was caused by the transient crack closure behavior shown 
schematically in Fig. 5, it was possible to calculate the crack closure stress 
intensity as a function of fatigue crack wake after the transition. The 
calculations were performed by first solving Eq. (3) for the values 

required to obtain the measured FCP rates after the transition to a higher 
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AK (Fig. 8). The obtained AK cc values were subtracted from K to 
obtain the K CL values. The results of the calculation are shown in Fig. 11 
and show equilibrium K^ is reached after 0.2 mm of crack growth. These 
results agree fairly well with the previously described surface removal test 
results which showed that can be fully reestablished after 0.125 mm of 

crack growth in case of plasticity induced closure (Fig. 4). 

The difference in the measured K values of the 3.2 mm the thick 

specimen between the lower and higher applied stress intensities is 

approximately 1 MPa/m as shown in Table 4. This difference in K values is 

considerably smaller in comparison to the thin specimens. This result is not 

very surprising considering that for the thick specimens both the lower and 

higher applied AK's resulted in a predominant plane strain state of stress. 

The predicted FCP rates Immediately after the transition from a low to higher 

AK as well as the predicted FCP rates after a plateau has been reached were 

calculated using Eq. (3). The predicted FCP rates were respectively 
-7 -7 

5.9x10 m/cycle and 4.6x10 m/cycle Immediately after transition and in 
the plateau region. The difference between these two values is only 30 percent 
as opposed to 130 percent for the thin specimen. 

No initial acceleration of FCP rates was detected after the transition to 


a higher AK even though a 30 percent Increase was predicted. However, it 


should be noted that the small scatter band present in obtaining the K 

w L 

values can have a significant effect on the predicted FCP rates. Also, the 
20 percent wide scatter band of the striatlon spacing measurements might have 
obscured the small Initial acceleration. Thus, since the gradient In 
after transition is small, it is possible that a small acceleration of FCP 
rates did take place after the transition to a higher AK, as predicted by the 
closure hypothesis, however, due to the scatter of the data no acceleration 
was detected. 
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CONCLUSIONS 


(1) Crack closure stress intensity was shown to increase at a 
significantly higher rate with respect to the crack wake length for plasticity 
Induced closure in comparison to roughness Induced closure. 

(2) For plasticity induced closure, no reduction in closure magnitude was 
found by removal of the crack wake to within 0.2 mm of the crack tip. For 
roughness induced closure, a reduction in closure magnitude was noted at a 
remaining crack wake length of 3 mm. 

(3) Roughness induced closure for specimen tested at AK of 2.5 MPa/m 
was found to be uniform throughout the thickness of the specimen, while near 
surface closure levels were approximately 50 percent higher than mid-thickness 
closure for specimens tested at AK of 8.8 and 17.6 MPayfn. 

(4) The 1.1 mm thick specimens exhibited a transient region of higher 
than expected FCP rates immediately after the transition from a low to high 
stress intensity. However, for the 3.2 mm thick specimen, no such load 
Interaction effects were detected. 

(5) This difference in transient behavior was explained by comparing K CL 
values for the lower and higher stress intensities. For the thin specimens, 
there was a large difference in K CL values between the lower and higher 
stress intensities due to change in state of stress resulting in a 
progressively decreasing AK g ^ as a function of crack length from transition. 
For the thick specimens, in a predominant plane strain state, the difference 

in K cl was small creating negligible variation in AK^ after transition 
to a higher aK. 
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TABLE 1 


TENSILE PROPERTIES - 7075-T6 


Orientation 

Ultimate 

strength 

0.2 percent 

Elongation, 


MPa 

ksi 

yield strength 

percent 




MPa 

ksi 


Longitudinal 

565 

82 

524 

76 

12 

Long transverse 

579 

84 

517 

75 

13 


TABLE 2 - TEST MATRIX 


Type of Test 

Thickness, 

mm 

Specimens 

Baseline data base - 

3 2 

4 

constant load amplitude 
Crack wake removal at: 



AX = 2.5 MPa/m 

3.2 

1 

AK = 8.8 MPa/m 

3 2 

2 

AK = 17.6 MPaTm 

3 2 

2 

Surface removal at: 



AK = 2.5 MPaTm 

3.2 ~ 1 .8 

1 

AK = 8 8 MPa/m 

3 2 - ~ 1.8 

1 

AK = 17 6 MPa/m 

3.2 -*-1.8 

1 

Low-high load sequence 

1.1 

2 

Low-high load sequence 

3 2 

2 

K ci at 



AK of 4 4 MPa Jm 

3.2 

1 

AK of 13.4 MPa-fm 

3 2 

1 

AK of 6 6 MPaw^tn 

1 1 

1 


AK of 16 MPa/m 





















TABLE 3. - PLASTIC ZONE SIZE 


AK, 

Plastic 

Plastic zone 

MPa/m 

zone size, 3 

size ratio 


mm 


2.5 

b 0.0087 

1 

8.8 

.11 

12.5 

17.6 

.438 

50 


Calculated by using r p = 1/* (K max /<> ys ) 2 
b Un1ty zone size 


TABLE 4 - MEASURE K QL VALUES 


Thickness , 

AK applied, 

‘a- 

mm 

MPa/m 

MPa/m 

3.2 

4.4 

1.54 

3.2 

13.4 

2.56 

1.1 

6.6 

2.3 

1.1 

16.0 

5.3 














THICKNESS, mm 

Figure 3. - K c , as a function of specimen 
thickness after surface removal. 



THICKNESS, mm 


Figure 4. - Increase in K C( _ after reinitiation 
of crack growth . 
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da/dn = 1x10-6 m/cycle 
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